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ABSTRACT 

The excess of antiprotons (P) observed in cosmic rays has been 
attributed to their production in super nova_ (SN) envelopes 
expanding in dense clouds. While creating P, y-rays are also 
produced and these clouds would shine as y-ray sources. The 
evolution of the y-ray spectrum has been calculated for clouds 
of n H = K)4 and lO^/cm 3 - 

1. Introduction. The observed flux of P in cosmic radiation [1-3] has 
invoked many new ideas in recent years to explain these observations [4] . 
The observed excess of P implies large amount of matter to be traversed 
by cosmic rays and recently, Mauger and Stephens [5] have suggested the 
possibility of producing P in the envelopes of SN exploding in dense 
clouds. Cosmic rays while traversing matter produce pions and kaons, in 
addition to P. These unstable particles decay to y-rays and electrons; 
y-rays are also produced by the interaction of electrons with matter and 
radiation fields. Such sources should shine in y-rays and it is neces- 
sary’ to look for them [6]. In this paper, we derive the evolution of 
Y-ray spectrum in SN envelopes, which expand in dense clouds, and examine 
the consequences. 

2 . Theoretical Approach . The evolution of nucleon and electron compo- 
nents in supernova envelopes can be examined by solving the following 
coupled equations 
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In the above equations, the 1st term on the R.H.S. describes the conti- 
nuous energy loss of particles. In the case of protons, this energy 
loss corresponds to ionization and adiabatic cooling, the latter being 
(dE/dt) A » l (E+2m)E/r (E+m) } (dr/dt) . The radius r and its derivatives 
are obtained from the dynamics of SN. The 2nd term describes the 
energy shift due to the finite elasticity of the interacting particle 
and the 3rd term corresponds to the loss of particles due to interaction. 
In these terms, X is the interaction mean free path and v the velocity 
of the interacting particle. In the case of electrons (Eqn.2), the 
continuous loss term contains also loss due to bremsstrahlung, inverse 
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Compton and synchrotron processes. For the inverse Compton process, we 
have used a radiation density corresponding to an optical outburst of 
10 ergs/s soon 1 after the SN explosion, which then decay with an e— 
folding time of 0.2 yi . The magnetic field inside the remnants is assu- 
med to scale as B 2 ot n H , with B = 4pG at n H = 1 atom. cm -3 . 

Equation 2 is coupled to the 1st equation through the term Q which 
is given by e 

O - r r , , dE <JE 3 

6 E y E tt e p 6 V/ ' {J P P v} dE p - ^ 2Tr(E Pj. d 9 

* • • (3) 

Here, rne integral over 9 describes the production of pions, in which p 
is the transverse momentum of pions, 9 the angle of emission and 
(E d 3 a/dp3) the invariant cross section, which depends upon E . The 
integrals over E^ and E^ take care of the energy distribution^of muons 
and electrons during decay; ¥ ' s are normalized functions. We have also 
included the knockon electrons, which are important below a few tens of 
MeV. All these parameters are taken from the work of Badhwar and 
Stephens [7] . The set of coupled equations (1) and (2) has been solved 
by Runge Kutt a method, it is assumed that the acceleration is complete 
at onset of adiabatic phase and the initial energy spectrum is assumed to 
be a power law in rigidity of the type A R~ 2 ’ 75 , where A = 2.5xl0 4 * /(m 2 . 
sr.s.GV/c) for nucleons and = 150/(m 2 .sr.s.GV/c) for electrons. The 
parameters relating to the evolution of SN in dense clouds have been 
described earlier [8], For the sake of simplicity, it is assumed that 
the size of cloudlets is such that the total amount of matter traversed 
by cosmic rays, when the envelope leaves the cloudlet, is about 50 g.cm -2 . 

During the expansion of the remnant, cosmic rays interact with 
matter to produce neutral pions which decay into Y-rays. Electrons 
interact to produce bremsstrahlung y-rays; the contribution from inverse 
Compton is very small during the adiabatic phase, y-rays thus produced 
are calculated using the following integrals. 

2 <3E .3 

p y (t) 7r° = 47T ; e f E / e ~T~ * {J D pv} d V {27T (E p, d 0> •••< 4 > 

TT° p Y * p dp 


P y (t) B = 4TT / e $ (E e ,E) py J e dE c ... (5) 

c 

Eqn. 4 is similar to Eqn. 3, except that tt° decays to y-rays and the 
energy distribution of y-rays is taken care through the integral over 
• The cross-section for tt° production is taken from Stephens and 
Badhwar [9] ; the bremsstrahlung cross-section 4>(E ,E) is with and with- 
out screening [10] . e 

3 • Y-ray Spectral Evolutxon. We have used for the calculation the in- 

terstellar cosmic ray spectrum to be the source spectrum in the SN, and 

to determine the total cosmic ray energy density in the source, the 

following procedure has been adopted. It is found that one needs 30% of 
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ENERGY IN GeV TIME IN YEARS 


Fig.l. y -ray spectra through tt° Fig. 2 y -ray brightness is plotted 

decay and bremsstrahlung are as a function of time for two n H 

shown for SN exploding in a values, 

cloud with 10 5 atom.cm~3 „ 

the cosmic ray nucleons to originate from SN in dense clouds in order to 
explain the P observations [8] . The remaining 70% is assumed to come 
from SN exploded in ordinary clouds with n^ * 10 atom.cm"^; the observed 
birth rate of these SN is one in ~30 yrs. Considering a galactic space 
with radius 15 kpc and thickness 0.5 kpc r one requires an energy release 
- 7xl0 60 eV in cosmic rays by a SN to account for the present energy 
density over a period of 3x10? yrs. It is assumed that the acceleration 
is complete in these sources around 200 yrs. The adiabatic cooling is 
expected to cease at the end of adiabatic phase, when the envelope frag- 
ments. Taking into account the energy loss processes during expansion/ 
the calculated energy output just after the acceleration is ^10^ 2 eV. 
Therefore, we use this factor to obtain y-ray brightness in our calcula- 
tions. 

The evolution of y-ray spectrum is shown in Fig. 1 for various 
stages of SN evolution in clouds of density n H = 10^ atom.cm"^. The 
ordinate brightness scale is given as photon per (GeV.s). It can be 
seen that during the early phase, the spectrum below a few hundred MeV 
is dominated by bremsstrahlung radiation and as a result, the total spec- 
trum can be represented by a simple power law. However, at the later 
stages, tt° decay y-rays become dominant. In order to examine the varia- 
tion of the total intensity with time, we have plotted in Fig. 2, the 
integral brightness above 100 MeV as well as above 30 MeV as a function 
of time, for n H = 10^ and 10^ atom cirT^; the upper scale is for n H = 10^ 
atom. cm" 3. 
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It can be noticed from Fig. 2 that, these sources are visible even 


upto far end of the Galaxy, if the threshold for detection is good to 
10” 6 photon/ cm 2 for n H = 10 3 atom. cm - 3 . The life time of this source is 
only a few hundred .yrs. and hence the number of such sources would be 
small: in a given time, in the case of n H = 10 4 atom. cm -3 , the threshold 
needed®detect them over the entire galaxy is <10“ 7 Photon. cm -2 . The 
energy dilution of cosmic rays in these sources from the onset of the 
adiabatic phase till they travel 50 g.cm -2 of matter is estimated to be 
by a factor 24.5 and 55 for njj = 10 3 and 10 4 atom. cm -3 respectively, 
which can be compared to 17.5 for n H = 10 atom. cm” 3 . Therefore, the 
number of such sources in the galaxy in comparison with the number of SN 
exploding in normal cloud is 0. 6:1.0 and 1.34:1.0 for n H = 10 3 and 10 4 
atom. cm -3 respectively. This would mean that equal number of SN are 
exploded in dense clouds with diameter $ 1 pc, but only a few of them 


would be visible for observation. 
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